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Abstract 

The production of neutral strange hadrons is investigated using deep-inelastic scattering 
events measured with the HI detector at HERA. The measurements are made in the phase 
space defined by the negative four-momentum transfer squared of the photon 2 < Q 2 < 
100 GeV 2 and the inelasticity 0.1 < y < 0.6. The K® and A(A) production cross sections 
and their ratios are determined. production is compared to the production of charged 
particles in the same region of phase space. The A — A asymmetry is also measured and 
found to be consistent with zero. Predictions of leading order Monte Carlo programs are 
compared to the data. 
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1 Introduction 



The production of strange hadrons in high energy particle collisions allows the investigation 
of strong interactions in the perturbative and non-perturbative regimes. Strange quarks are 
created in the non-perturbative process of colour string fragmentation, which constitutes the 
dominant production mechanism of strange hadrons. In deep-inelastic scattering (DIS), strange 
quarks also originate from the strange sea in the nucleon, boson-gluon fusion and heavy quark 
decays. Measurements of strangeness production have been used to investigate the suppression 
of strangeness relative to lighter flavours in fragmentation. The universality of fragmentation 
in different processes can be studied by comparing differential cross sections of the production 
of K® and A (A) hadrons in various regions of phase space. Further information is gained by 
studying the ratios of production rates of A (A) to K® and of K° to charged hadrons (/i ± ) as 
some model dependencies are expected to cancel. 

The baryon production mechanism was studied in e + e~ annihilation [[Ill^ElHEl, a process 
without incident baryons. Data involving a baryon in the initial state, like ep collisions at 
HERA, provide additional information. In particular, data on the A - A production asymmetry 
from HERA are of interest as an experimental constraint for theories of baryon number transfer 
[0. Fixed target data have shown [7] that the A production rate substantially exceeds that of the 
A in the so-called remnant region because the baryon number of the target is conserved. 

This paper presents a measurement of neutral strange particle (K® and A) production in DIS 
at negative four momentum transfer squared 2 < Q 2 < 100 GeV 2 and at low values of Bjorken 
x. The study is based on data collected with the HI detector at HERA at a centre-of-mass 
energy of 319 GeV in the years 1999 and 2000. This data sample is 40 times larger than that 
used in the previous HI publication [8 J and covers a wider kinematic range. Measurements 
of and A production in different kinematic ranges have also been reported by the ZEUS 
collaboration [f9j. The differential cross sections of K° mesons, A (A) baryons and their ratio 
as well as the ratio of K® to charged hadrons are presented as a function of various kinematic 
variables, both in the laboratory frame and in the Breit frame. The results are compared with 
predictions obtained from leading order Monte Carlo calculations, based on matrix elements, 
with parton shower simulations. The main feature of the data is a suppression of strange quark 
production relative to lighter quarks; this is discussed within the context of the framework of 
the LUND El fragmentation model. 



2 Phenomenology 

2.1 Production of Strange Hadrons 

Particles with strangeness can be produced in DIS in the hard sub-process and in the hadronisa- 
tion of the colour field, as illustrated schematically in figured] 

Figure \T^l) shows strangeness production within the quark parton model (QPM), where a 
strange quark s from the nucleon sea participates in the hard interaction. Figured}?) illustrates 
s production in a boson-gluon fusion (BGF) process, where a gluon emitted from the nucleon 
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Figure 1: Schematic diagrams of the different processes contributing to strangeness production: 
a) direct production in the QPM, b) BGF, c) decays of heavy quarks and d) hadronisation. 



splits into a ss quark pair. Figure [Tfc) depicts heavy quark (charm c and beauty b) production 
by boson-gluon fusion (BGF) with subsequent weak decay into s(s) quarks. This process is 
suppressed at low Q 2 due to the masses of the heavy quarks. These production mechanisms 
(figures b, c) are characterised by a hard scale allowing for a perturbative treatment. The 
relative rate of the BGF processes depends strongly on the Bjorken scaling variable x due to the 
strong rise of the gluon density at low x. In the kinematic region studied in this analysis (low x) 
the BGF contributions are expected to be significant. According to the Monte Carlo predictions 
described below, roughly 25 % of the strange hadrons originate from strange quarks produced 
in the hard interaction either directly (figures [T^ and b) or through heavy quark production in 
BGF processes (figure Ht). In regions of phase space where the quark masses are not relevant 
with respect to the process scales (e.g. at very high Q 2 ) this rate can reach up to 50%. 

The largest contribution to strange quark production is due to the colour field fragmentation 
processes, as illustrated in figure[0i). As these processes occur at large distances they cannot be 
treated perturbatively and thus phenomenological models, such as the LUND string model ifTOll . 
are required for their description. 

Frames of reference customarily used to study particle production are the laboratory and the 
Breit frame [fTTTl . In the Breit frame of reference the virtual space-like photon has momentum Q 
but no energy. The photon direction defines the negative ,2-axis with the proton moving in the +z 
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direction. The transverse momentum in the Breit frame pff rett is computed with respect to this 
axis. Particles from the proton remnant are almost collinear to the incoming proton direction, 
therefore the hemisphere defined by p^ reit > is labelled as the target hemisphere. Equally, 
in the QPM the struck quark only populates the current hemisphere (p^ reit < 0). Higher order 
processes modify this simple picture as they generate transverse momentum in the final state 
and may lead to particles from the hard subprocess propagating into the target hemisphere. 

In the current hemisphere, the mechanism of particle production should in principle resem- 
ble that of collisions without an incident proton like e + e~. In analogy with e + e~ collisions the 
fragmentation variable x® reit = 2\p \/Q is defined, where p is the momentum of the particle 
in the Breit frame; x^ reit corresponds to x p = p/pbeam in e + e~ collider experiments. Strange 
quarks produced directly in the hard interaction are expected to preferentially populate the cur- 
rent hemisphere, which is less sensitive to non-perturbative strangeness contributions. In the 
case of baryon production the hemisphere separation is useful to study also baryon transfer, 
which is expected to be relevant at high x^ reit in the target frame. 



2.2 Monte Carlo Simulation 

The deep-inelastic ep interactions are simulated using the DJANGOH program 03 . It gener- 
ates hard partonic processes at Born level and at leading order in as (e.g. 7*g — > q, 7*g — > qg, 
1*9 ~^ QQ X convoluted with the parton distribution function (PDF) for the proton, chosen herein 
to be CTEQ6L |fT3ll . The factorisation and renormalisation scales are set to fij = = Q 2 . 
Within DJANGOH, higher order QCD effects are accounted for using either the parton shower 
approach as implemented in LEPTO 03J (referred to as MEPS) or by the so-called colour 
dipole model approach available within ARIADNE O (referred to as CDM (161). In LEPTO, 
the parton showers are ordered in the transverse momenta (kx) of emissions, according to the 
leading \og(Q 2 ) approximation. In the ARIADNE program, the partons are generated by colour 
dipoles spanned between the partons in the cascade; since the dipoles radiate independently, 
there is no k T ordering. 

The hadronisation process is modelled according to the LUND colour string fragmentation 
model ifTOl . as implemented in the JETSET ifTTl program. Within this model, the strange quark 
suppression is predominantly described by the (constant) factor X s = P s /P q , where P s and P q 
are the probabilities for creating strange (s) or light (q = u or d) quarks in a non-perturbative 
process from the colour field during the fragmentation process. Further important parameters of 
this model are the diquark suppression factor X qq = P qq /P q , i.e. the probability of producing a 
light diquark pair qqqq from the vacuum with respect to a light qq pair, and the strange diquark 
suppression factor \ sq = (P sq /P qq )/(P S / P q ), which models the relative production of strange 
diquark pairs. These are the two most relevant factors for the description of baryon production. 
The ss pair production rate is primarily dominated by A s , i.e. u(u) : d(d) : s(s) — 1 : 1 : A s . 
The values tuned to hadron production measurements by the ALEPH collaboration [4] (\ s = 
0.286, X qq = 0.108, and X sq = 0.690) are taken herein as default values for the simulation of 
hadronisation within JETSET. 

Previously published HI and ZEUS data [IHEEBl are better described by a lower value X s = 
0.2. A recent ZEUS analysis [9] favours X s = 0.3 from cross section results and X s = 0.22 from 
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measurements of the strange mesons to charged hadrons ratio. The same theoretical framework 
is also used in e + e~ analyses and thus allows for tests of strangeness suppression universality. 

Monte Carlo event samples generated with DJANGOH are used for the acceptance and 
efficiency correction of the data. All generated events are passed through the full GEANT 
[fT9l based simulation of the HI apparatus and are reconstructed and analysed using the same 
programs as for the data. 



3 Experimental Procedure 

3.1 HI Detector 

A detailed description of the HI detector can be found in EOl . In the following, only those 
detector components important for the present analysis are described. HI uses a right handed 
Cartesian coordinate system with the origin at the nominal ep interaction point. The proton 
beam direction defines the positive z-axis of the laboratory frame and transverse momenta are 
measured in the x — y plane. The polar angle 9 is measured with respect to this axis and the 
pseudorapidity 77 is given by i] = — In tan |. 

Charged particles are measured in the Central Tracking Detector (CTD) in the range — 1 . 75 < 
i] < 1.75. The CTD comprises two cylindrical Central Jet Chambers (CJCs), arranged con- 
centrically around the beam-line, complemented by a silicon vertex detector (CST) [21 J, two 
,2-drift chambers and two multi-wire proportional chambers for triggering purposes, all within a 
solenoidal magnetic field of strength 1.16 T. The transverse momentum resolution is a(px) /pt 
~ 0.006 p T I GeV © 0.015 E2l . In each event the tracks are used in a common fit procedure to 
determine the ep interaction vertex. 

The tracking detectors are surrounded by a Liquid Argon calorimeter (LAr) in the forward 
and central region (—1.5 < r] < 3.4) and by a lead- scintillating fibre calorimeter (SpaCal) 
in the backward region [23] (—4 < i] < —1.4). The SpaCal is designed for the detection 
of scattered positrons in the DIS kinematic range considered here and has an electromagnetic 
energy resolution of a E /E ~ 7%/^/E/ GeV © 1%. The backward drift chamber (BDC), 
positioned in front of the SpaCal, improves the measurement of the positron polar angle and is 
used to reject neutral particle background. The DIS events studied in this paper are triggered by 
an energy deposition in the SpaCal, complemented by signals in the CJCs and in the multi-wire 
proportional chambers. 

The luminosity is determined from the rate of the Bethe-Heitler process, ep — > epj, mea- 
sured using a calorimeter located close to the beam pipe at z = —103 m. 

3.2 Selection of DIS Events 

The analysis is based on a data sample corresponding to an integrated luminosity of C = 
49.9 pb _1 , recorded when HERA collided positrons at an energy E e = 27.6 GeV with protons at 
920 GeV in the years 1999 and 2000. 
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The selection of DIS events is based on the identification of the scattered positron as a 
compact calorimetric deposit in the SpaCal. The cluster radius is required to be less than 3.5 cm, 
consistent with an electromagnetic energy deposition. The cluster centre must be geometrically 
associated with a charged track candidate in the BDC. These conditions reduce background 
from photoproduction processes. 

At fixed centre of mass energies y/s the kinematics of the scattering process are described 
using the Lorentz invariant variables Q 2 , y and x. These variables can be expressed as a function 
of the scattered positron energy E' e and its scattering angle 6 e in the laboratory frame: 

g» = 4** »»(|), » = i-|w(!), , = g. (i) 

The negative four-momentum transfer squared Q 2 and the inelasticity y are required to lie 
in the ranges 2 < Q 2 < 100 GeV 2 and 0.1 < y < 0.6. Background from events at low Q 2 , 
in which the electron escapes undetected down the beam pipe and a hadron fakes the electron 
signature, is suppressed by the requirement that the difference H(E — p z ) between the total 
energy and the longitudinal momentum must be in the range 35 < Y>(E — p z ) < 70 GeV, where 
the sum includes all measured hadronic final state particles and the scattered electron candidate. 
Events are accepted if the 2-coordinate of the event vertex, reconstructed using the tracking 
detectors, lies within ±35 cm of the mean position for ep interactions. 

3.3 Selection of Hadron Candidates 

The neutral strange K® meson and A baryon states^] are measured by the kinematic reconstruc- 
tion of their decays K® — > n + ix~ and A — > pn~ . The analysis is based on charged parti- 
cles measured in the central region of the HI detector with a minimum transverse momentum 
Pt > 0.12 GeV. The neutral strange hadrons K® and A are identified by fitting pairs of oppo- 
sitely charged tracks in the x — y plane to their secondary decay vertices, with the direction of 
flight of the mother particle constrained to the primary event vertex. and A candidates are 
retained if the fit probability is above 1 %. In order to reduce background, the radial distance 
L of the secondary vertex to the beam line is required to be larger than 5 mm and the vertex 
separation significance L/gl > 4, where gl is the uncertainty of L. The transverse momentum 
and the pseudorapidity of the K® (A) candidates are required to satisfy 0.5 < pr < 3.5 GeV 
and \t)\ < 1.3. A detailed description of the analyses can be found in If24ll25ll . 

For K° candidate reconstruction both tracks are assumed to be pions, while for the A re- 
construction the track with the higher momentum is assumed to be the proton and the other 
track is assumed to be the pion. The contamination from A (K®) decays in K®(A) candidates is 
suppressed by a rejection of the corresponding invariant mass region: | M(irp)— |> 6MeV 
for the K® and | M^ttti) — m K a |> 10 MeV for the A selection. The A (A) baryons are tagged 
by the electrical charge of the decay proton (antiproton). The invariant mass spectra M(n + n~) 
and M(pix) of all candidates passing these criteria are shown in figures [2] and [3l respectively. 

'Unless explicitly mentioned, a reference to a state implicitly includes the charge conjugate of that state. 
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The number of signal particles N s is obtained by fitting the invariant mass spectra with the 
sum of a signal and a background function. The signal function S has the same shape for K® and 
A and is composed of two Gaussian functions of identical central value \i and of different widths 
<Ti and o"2 that account for different resolution effects. The background functions B K o(M) and 
B\(M) are chosen with different shapes for the K° s and A cases. These functions are defined 
according to 

S(M) = P -G(Ns,fi,a 1 ) + (l-P )-G(N s ,fi,a 2 ), (2) 
B K o(M) = P, + P 2 -M, (3) 
B A (M) = P x -(M-m K ) p > . e (i+*VM+P4.M3)_ (4) 



Here, M denotes the nn and the pn invariant mass, respectively, and rriA the nominal mass 
of the A baryon ll26ll . The normalisation Ns, the central value /x, the widths a\ and a 2 of the 
Gaussian function G and the parameters Pj are left free in the fit. P represents the relative 
normalisation of the two signal Gaussians. For the differential distributions the fit is repeated in 
each of the kinematic bins. 

The fit yields approximately 213000 K° s mesons. The fitted mass of 496.9 ± 0.1 ( stat.) MeV 
is consistent with the world average Il26l and the measured mean width 13.8 ± 0.4 ( stat.) MeV 
is described by the simulated detector resolution within 20%. In the case of the A the fit yields 
approximately 22000 A and 20000 A baryons. The fitted mass of 1115.8 ± 0.1 ( stat.) MeV is 
also consistent with the world average ll2~6ll and the measured width of 4.3 ± 0.3 ( stat.) MeV 
is consistent with the detector resolution within 20%. 

Charged hadrons /i ± used for the ratio R^K^/h^ are defined as long-lived particles with 
a lifetime > 10~ 8 s detected in the same kinematic region as strange particles (\t]\ < 1.3, 
0.5 < pr < 3.5 GeV), with the following additional requirements: each track must point to 
the primary vertex, the number of associated hits in the CJC must be greater than eight, the 
radial track length must be longer than 10 cm and the radial distance from the beam line to the 
innermost hit associated with the track must be less than 50 cm. 



4 Results 



4.1 Determination of Cross Sections 

The total inclusive cross section a V i S in the accessible kinematic region is given by the following 
expression: 

a ms (ep^e[K° s ,A,h ± }X) = - / , (5) 

C ■ e ■ BR ■ (1 + d rad ) 

where N represents the observed number of K®, the sum of A and A baryons or the charged 
hadrons h ± , respectively. C denotes the integrated luminosity. The branching ratios BR for the 
and A decays are taken from [26J and BR = 1 for charged hadrons. The number of and 
A particles are determined by fitting the mass distributions as explained in section 13.31 In the 
case of differential distributions the same formula is applied in each bin. 
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The efficiency e is given by e = e rec -e tr ig, where e rec is the reconstruction efficiency and e trig 
is the trigger efficiency. The reconstruction efficiency is estimated using CDM Monte Carlo 
event samples for the kinematic region and the visible range defined in sections 13.21 and 13 .31 
and amounts to 33.3 % and 19.5 % for the mesons and the A baryons, respectively. These 
numbers include the geometric acceptance and the efficiency for track and secondary vertex 
reconstruction. The geometric acceptance to find both decay particles in the CTD is about 80 % 
for the mesons and 70 % for the A baryons, respectively 

The trigger efficiency is extracted from the data using monitor triggers and amounts to 
81.5 % and 83.3 % for the K° s and the A, respectively. The radiative correction 5 ra( i corrects the 
measured cross section to the Born level and is calculated using the program HERACLES ll27ll . 
It amounts to 5 ra d = 6.6(4.3) % for the K® (A) on average and varies between —8% and +19% 
over the kinematic range considered. The trigger efficiency and radiative corrections are as- 
sumed to be the same for particles and antiparticles. 

In the case of charged hadrons h ± , the reconstruction efficiency e rec is defined such that it 
includes corrections for K° and A decays, secondary interactions, photon conversions and the 
track reconstruction efficiency. The total correction e(l + S ra d) amounts to 81.1%. 

4.2 Systematic Uncertainties 

The systematic uncertainties are studied using the CDM Monte Carlo simulation, unless oth- 
erwise stated. For the inclusive cross sections, the resulting systematic uncertainties are sum- 
marised in tabled! For the differential cross sections, the systematic uncertainties are estimated 
in each bin. The following contributions are considered: 

• The energy scale in the Spacal measurements is known to 1 %, except for the lowest Q 2 
bin (2 < Q 2 < 2.5 GeV 2 ) where the uncertainty on the energy measurement is 2.5 %. 

• The uncertainty of the measurement of the polar angle of the scattered positron is 1 mrad. 

• The uncertainty on the overall number of reconstructed K® and A particles is determined 
from data by comparing the numbers obtained from the fit of the mass spectra with the 
number obtained by simply counting the events within ±6 a of the nominal mass after 
subtracting the expected background. The number of background events is estimated 
by integrating the background function described in equation |4] over the corresponding 
interval (0.42 — 0.58 GeV for K® and 1.085 — 1.2 for A). The procedure is cross checked 
by performing the fit in different mass ranges. 

• The uncertainty of the reconstruction efficiency is determined by comparing its estimation 
using different models. The uncertainty is taken as 50 % of the difference between the 
CDM and the MEPS Monte Carlo simulations. 

• The uncertainty on the trigger efficiency is obtained by comparing estimates using differ- 
ent monitor triggers (MT). 

• The uncertainty on the luminosity measurement is 1.5 %. 
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• The uncertainty of the charged hadron reconstruction is 2 % per track. For the measure- 
ment of the A to K® ratio the uncertainty caused by the pion track appearing in both 
decays is assumed to cancel. The systematic uncertainty on the ratio / is estimated 
to be 2.0 %. 

• The uncertainty due to the decay branching ratios is taken as 0.8% for A and is negligible 
for K° s m. 
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±0.4 


±1.2 


±1.2 


±3.5 


trigger efficiency 
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trig trig 

MTsetl 

trig 


+0.4 
-0.9 


+1.0 

-1.4 


+1.1 

-1.6 


+0.4 
-1.0 


luminosity 




±1.5 


±1.5 






track reco. 


2.0% per track 


±4.0 


±4.0 


±2.0 


±2.0 


branching ratio 




±0.1 


±0.8 


±0.8 


±0.1 


Total systematic uncertainty 


+5.6 
-5.8 


+5.8 
-6.0 


+3.1 

-3.3 


+4.1 
-4.2 



Table 1: Systematic sources, variations and corresponding relative errors of the inclusive cross 
sections and of the ratios of A to and to charged hadrons. All relative errors are expressed 
as percentages. 



The systematic errors due to these uncertainties are estimated by varying each quantity 
within its error in the Monte Carlo simulation and repeating the cross section measurement. 
In the cross section calculation, the contributions are added in quadrature and included in the 
uncertainty shown in the individual bins of the differential distributions. In the ratios, the un- 
certainties on the electron energy scale and polar angle, as well as the luminosity, cancel. The 
other sources of uncertainty are assumed to be uncorrelated and are added in quadrature. 

4.3 Inclusive Production Measurements 

The inclusive K®, A and charged hadron hr^ production cross sections a vis are measured in the 
kinematic region 2 < Q 2 < 100 GeV 2 and 0.1 < y < 0.6, for the ranges 0.5 < Pt(K q s} A, /i ± ) < 
3.5 GeV and \r](K®,A, /i ± )| < 1.3. 
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The K® cross section is found to be 

Vvis{ep -> eK° s X) = 21.18 ± 0.09(stat)^;51(syst) nb. (6) 

The measurement is in agreement with the expectation 21.77 nb based on the LO Monte Carlo 
program DJANGOH, using the CDM approach and the default value of X s = 0.286. 

The cross section for the sum of A and A baryon production is measured in the same kine- 
matic region and is found to be 

cr vis {ep -> e[A + A]X) = 7.88 ± 0.10(stat.)+°;^(syst.) nb, (7) 

in agreement with the expectation of 7.94 nb from the DJANGOH calculation. The individual 
A and A production rates are measured to be 

Ovis{ep -> eAA) = 3.96 ± 0.06(stat.)+°;|[(syst.) nb, (8) 
a ms (ep eAX) = 3.94 ± 0.07(stat.)i°;|j(syst.) nb, (9) 

and are therefore found to be consistent with each other within the statistical precision. The 
measurements are also in agreement with the DJANGOH prediction of 3.97 nb. The systematic 
errors are fully correlated. 



The inclusive ratio of strange baryon to meson production is determined to be 

0.372 ± 0.005(stat.)l°;°^(syst.), (10) 



crvis(ep -> e[A + A]X) _ n _ , n (Wir ,,„„, +0 . riJ , , 



o- V is(ep -> eK° s X) 

in agreement with the prediction of 0.365 from the DJANGOH calculation. 
The ratio of cross sections of mesons to charged hadrons ^ is found to be 

cr vis {ep ->• eK®X) _ nnnoCotat ^ +0.0019, 



o- V is(ep -> eh±X) 



0.0645 ± 0.0002(stat.)+^^(syst.), (11) 



in agreement with the DJANGOH prediction of 0.0638 based on MEPS with A s = 0.22. 
Similar values of 0.05 — 0.07 are obtained for the ratio of the average A° multiplicity over the 
average charged pion multiplicity in e + e~ annihilation events at centre of mass energies from 
10 to 200 GeV [|26l. 



4.4 Differential Production Cross Sections 



Production cross sections and ratios of K®, A and charged hadrons are measured in the 
visible kinematic region differentially in the event variables Q 2 and x and in the laboratory 
frame variables px and 77. Differential cross sections are also measured as a function of the 
variables x^ reit and p^ reit defined in the Breit frame. The results are bin-averaged and no bin- 
centre corrections are applied. The distributions are shown in figures |4] to [T3] and are compared 
with the predictions. The numerical values are also listed in tables |2] to [U 
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4.4.1 Discussion of K° and A Results 



The measured differential cross sections of K® and A production are shown in figures |4]to|7] and 
listed in tables [2] to [5] The cross sections decrease rapidly as a function of Q 2 and x, similarly 
to the inclusive DIS distributions. The cross sections are also observed to fall rapidly with pt. 

In the laboratory frame the overall features of the distributions are reproduced by the DJAN- 
GOH simulations at the level of 10 to 20%. For comparison, the CDM and MEPS model pre- 
dictions are each given with two values of the suppression factor X s = 0.3 and \ s = 0.22. The 
predictions based on the CDM model with \ s = 0.3 provide a reasonably good description of the 
data for K® and A production. The MEPS simulation produces distributions, which are quite 
similar in shape to the CDM model predictions but with a different normalisation in the case 
of K® production, where a lower value of \ s = 0.22 describes the data better. In the case of A 
production, both MEPS and CDM predictions are very similar in shape and normalisation and 
X s = 0.3 provides a better description of the data. For these comparisons, only the parameter X s 
is varied to describe the data. However, in contrast to the K®, the A production cross sections 
also depend significantly on the JETSET parameters that describe diquark and strange diquark 
creation. 

The cross sections measured as a function of xl 3 reit and p^ Teit in the Breit frame are shown in 
figures [6] and [7] and listed in tables 0] and [51 for both the target and the current region. The cross 
section values in the target regions are about one order of magnitude higher than in the current 
region. They are generally well described by both the MEPS and CDM model predictions. The 
predicted momentum distributions tend to be softer than in the data. However, in the current 
region the sensitivity to \ s is clearly reduced with respect to the laboratory frame or the target 
region. This is due to both larger errors and an increased fraction of strangeness produced in 
perturbative processes, which contributes up to about 50% (compared to about 25% in the target 
hemisphere). 

To test the mechanism of baryon number transfer, the asymmetry in the production of A 
with respect to A is measured by the variable 

_ cr vis (ep -> eAX) - Qjnsjep -> eAX) 
Ovis{ep -> eAX) + a vis (ep -> eAX) ' 

A significant A - A asymmetry ^ would indicate a transfer of the baryon number from 
the proton beam to the final state strange particles. The measured distributions of in the 
laboratory and Breit frames are shown in figures [8] and [9l respectively. All distributions are 
observed to be compatible with zero within errors. Thus, no evidence of baryon number transfer 
is visible in the measured A/A data. 

In order to test for possible dependencies of strange hadron production on the proton parton 
density functions, the measured distributions are compared with different PDF parametrisa- 
tions. Figure [10] shows the differential cross sections for K® and A production compared with 
the CDM predictions using the CTEQ6L flU, H12000LO fll! and GRV LO flU parametrisa- 
tions and \ s = 0.286. The predictions of the Q 2 dependence of the cross section are notably 
different for different PDFs for both the and the A. The p T distributions indicate only a 
slight dependence while the 77 distributions do not exhibit any PDF dependence. The small 
discrepancy in the forward direction is not resolved by different PDF parametrisations. Similar 
results are obtained in the Breit frame. 
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4.4.2 Ratios of Production Cross Sections 

Different aspects of baryon production within the fragmentation models can be tested with 
reduced theoretical uncertainties by studying the ratio of the differential cross sections for A 
baryons and mesons R(A/ K®) = da(ep — > eAX) / da(ep — > eK®X). The measurements 
are shown in figures [TT] and [12] and listed in tables [6] and [71 The CDM implementation provides 
a reasonably good description of the data in the laboratory frame (figure [TTI). although system- 
atic deviations are seen at high Q 2 and in the shape of the r\ distribution, whereas the MEPS 
predictions clearly underestimate the data. The model predictions are not sensitive to X s , as 
expected. 

The dependence of R(A/K®) on p® re%t and x? reit (figure fT2l) are reasonably well described 
in both the target and current hemispheres. The predictions are almost independent of the model 
implementation (CDM, MEPS) and the A s values used. 

The ratio of differential production cross sections for mesons and charged hadrons, 
denoted by R^K^/h^) = da(ep — > eK®X) I da(ep — > efi^X), is equivalent to the ratio of the 
average multiplicities of and charged hadrons. In contrast to inclusive production, the 
correlation of R(K® //i ± ) to the parameter A s is expected to be less model dependent. By taking 
the ratio, inadequacies of the model description of the partonic final states and in particular 
of the dependence on the proton structure function should cancel to a large extent. The ratio 
R^K^/h^ is shown in figure[[3]and listed in table[8]as a function of Q 2 , x, pt and 77. The ratio 
strongly rises with increasing px and remains approximately constant as a function of all the 
other variables. This px dependence of R(K® / h^) reflects a general kinematic feature (heavier 
particles receive the larger fraction of the system momentum) and can also be observed in the 
AIK® ratio in figure dU 

Also shown with the data are the CDM and MEPS model predictions for two values of X s 
(0.22 and 0.3). Overall, no single prediction is able to fully describe the shapes of all R^K^/h^) 
distributions, failing in particular in the low p T , low x and large positive 77 regions. The pr 
spectrum of R( K® /h^) is found to be harder in the data, consistent with the conclusions derived 
from the cross section measurement. 

The shapes of the ratios RijC^/h^) are reasonably described by both CDM and MEPS 
model predictions. However, there is a difference in normalisation between the two models. 
The CDM prediction with A s = 0.3 is in better agreement with the data at low Q 2 , whereas at 
high Q 2 a value of A s = 0.22 is preferred, as observed in the ZEUS data [9|. In contrast, the 
MEPS model predictions prefer a lower value of A s = 0.22 over the full phase space. 

A comparison of the predictions, applying different settings for the diquark-quark suppres- 
sion factors (X qq , \ sq ) shows the expected behaviour. In general, no changes are visible in the 
shapes of the differential distributions, however some differences are present in the absolute 
normalisation. The K® distributions are not affected, as expected, and both the A and the ratio 
R(A/K®) show the anticipated correlations to the suppression factors (\ qq , X sq ). These pre- 
dicted effects are mostly independent of the choice of the A s value, used for the simulation, and 
indicate that the "ALEPH-tune" from e + e~ collisions also describes the overall features of the 
data in ep collisions, supporting the universality of strangeness production. 
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5 Conclusions 



The production cross sections and ratios of the production of K®, A and charged hadrons h 
are measured inclusively and also differentially as a function of the DIS variables and of the 
final state particle variables in the visible kinematic region, defined by 2 < Q 2 < 100 GeV 2 , 
0.1 < y < 0.6,0.5 <p T (K s °,A,/i ± ) < 3.5 GeV and \rf(K%, A, /i ± )| < 1.3. 

The measured total cross sections and their ratios are in agreement with the predictions 
based on DJANGOH. The overall features of the various differential distributions are reasonably 
well reproduced by both simulations, based either on the CDM or the MEPS approach, when 
applying model parameters obtained from e + e data at LEP. However, predictions based on 
a single value of \ s fail to describe the details of the distributions in various regions of the 
phase space, in particular in the low pr, low x and large positive 77 regions. The production of 

and A particles, as measured in the Breit frame, is in general described by both CDM and 
MEPS predictions. The measurement of the asymmetry in the production of A with respect to 
A, which is found to be consistent with zero within errors, does not support the hypothesis of 
baryon number transfer. 

The A to Kg cross section ratio is better described by the CDM prediction and is nearly 
independent of the A s value, whereas for the K® to charged hadrons cross section ratio the 
MEPS model with A s = 0.22 is in better agreement with the data. 
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Figure 2: The invariant mass spectrum for it + it particle combinations. The data are shown 
with error bars denoting the statistical uncertainty. 
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Figure 3: The invariant mass spectra for a) A — > pn and b) A — > p7r + particle combinations. 
The data are shown with error bars denoting the statistical uncertainty. 



18 



£ 7 
O 

I 6 



O 

73 



o w 

CD 

t 
CL 



H1 



5 a 



CO 
Q 

o 



a) 



ep 



H1 Data 
CDM (X s =0.3) 
CDM (X s =0.22) 
MEPS {X s =0.3) 
MEPS (X s =0.22) 




1.2 _ri_i ' 

0.8 ' 



10 



10 2 
Q 2 [GeV 2 ] 



eK°X 



f 10 = 



O (0 

T 

Q. 

TJ 



H1 



10 



4 - 



-f- H1 Data 

CDM (X s =0.3) 

CDM (>. s =0.22) 

MEPS (X s =0.3) 

----MEPS(X S =0.22) 



10 J 



b) 



S 1.2 



1=*= 



g 0.8 



=i=j ».., 



10' 



10" ; 



10- 

x 



> 

52 40 

i- 

Q. 

>< 

o « 

t 20 
a. 

73 



S 1.5 

(8 

a 



o 

CD 



1 



H1 



c) 



-f - H1 Data 

CDM (X s =0.3) 

CDM (>. s =0.22) 

MEPS (X s =0.3) 

----MEPS (X s =0.22) 




0.5 



1.5 



2.5 3 
P T [GeV] 



3.5 



O (0 

CD 

t 
Q. 
CD 

¥ 



10 

8 
6 
4 
2 




H1 Data 
CDM (X s =0.3) 
CDM (X s =0.22) 
MEPS (X s =0.3) 
MEPS (X s =0.22) 




11 



Figure 4: The differential production cross sections for K® in the laboratory frame as a func- 
tion of the a) photon virtuality squared Q 2 , b) Bjorken scaling variable x, c) transverse 
momentum p T and d) K® pseudorapidity r]. The inner (outer) error bars show the statistical 
(total) errors. On the bottom of each figure, the "Theory/Data" ratios are shown for different 
LO Monte Carlo predictions (see text). For comparison, the data points are put to one and only 
uncorrelated errors are shown; the correlated systematic errors are indicated by the grey band. 
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Figure 5: The differential production cross sections for A in the laboratory frame as a function 
of the a) photon virtuality squared Q 2 , b) Bjorken scaling variable x, c) A transverse momentum 
p T and d) A pseudorapidity rj. More details in the caption of figure HI 
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Figure 6: The differential production cross sections in the Breit frame as a function of 
transverse momentum p^, reit and momentum fraction xL 3 rett in the target hemisphere (a, b) and 
in the current hemisphere (c, d). More details in the caption of figure HI 



21 



> 
o 
O 
B 
c 



x 

< 

T 

Q. 

O 

¥ 



8 



H1 



ft 



a) 



« 1.5 



(0 
Q 

O 




1 • ±± 



e p — > e A X (Breit frame) 

hT 



H1 Data 
CDM (X s =0.3) 
CDM (X s =0.22) 
MEPS (?i s =0.3) 
MEPS (X s =0.22) 



target region 



0.5 1 1.5 2 2.5 3 3.5 4 

.Breit 



p Bre.t [QeV] 



H1 Data 

CDM (X s =0.3) 

CDM (>. s =0.22) 
- MEPS (X s =0.3) 
----MEPS(X S =0.22) 




,Breit 



> 
O 

O 

-Q 
C 



0.3 



x 

< 

CD 

T 

Q. 

To 

TJ 



(0 
Q 

O 




0.2 



0.1 



1.2 
1 

0.8 



H1 



c) 



0.5 



-f - H1 Data 

CDM (X s =0.3) 

CDM (>. s =0.22) 

MEPS (X s =0.3) 

----MEPS (X s =0.22) 



current region 



1.5 



2 2.5 3 
pBreit [GeV] 



CQ Q 

X 

x 

< 

T 

Q. 

"b" 

TJ 



CO 
Q 



O 

CD 



0.6 



0.4 



0.2 



1.2 
1 

0.8 



H1 

— n — 


-+- H1 Data 

CDM (X s =0.3) 

CDM (/. s =0.22) 




MEPS (X s =0.3) 

----MEPS(X S =0.22) 


d) 




f 

current region 








— 1 1 — ' 









0.2 0.4 0.6 0.8 



, Breit 



Figure 7: The differential production cross sections for the A baryons measured in the Breit 
frame as a function of A transverse momentum p^ reit and momentum fraction xi 3 re%t in the 
target hemisphere (a, b) and in the current hemisphere (c, d). More details in the caption of 
figure HI 
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Figure 8: The asymmetry A a of the differential production cross sections of the A and A baryons 
in the laboratory frame as a function of the a) photon virtuality squared Q 2 , b) Bjorken scaling 
variable x, c) transverse momentum p T and d) pseudorapidity rj. The asymmetry is defined as 

A A = [cr vis (ep -> eAX) - a vis (ep -> eAX)]/[cr™ s (ep -> eAX) + a vis (ep -> eAX)]. The error 
bars show the statistical uncertainty. 
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Figure 9: The asymmetry A\ of the differential production cross sections of the A and A baryons 
measured in the Breit frame as a function of transverse momentum p^, rett and momentum frac- 
tion Xp reit in the target hemisphere (a, b) and in the current hemisphere (c, d). The asymmetry is 
defined as A A = [cr vis (ep — > eAX) — a vis (ep — > eAX)]/[a vis (ep — > eAX) + a vis (ep — > eAX)]. 
The error bars show the statistical uncertainty. 
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Figure 11: The ratio of the differential production cross sections for A baryons and K° mesons 
in the laboratory frame as a function of the a) photon virtuality squared Q 2 , b) Bjorken scaling 
variable x, c) transverse momentum px and d) pseudorapidity 77. More details in the caption of 
figure HI 
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Figure 12: The ratio of the differential production cross sections for A baryons and mesons 
in the Breit frame as a function of transverse momentum p^ reit and momentum fraction xp 3 reit 
in the target hemisphere (a, b) and in the current hemisphere (c, d). More details in the caption 
of figure HI 
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Figure 13: The ratio of the differential production cross sections for K® mesons and charged 
hadrons in the laboratory frame as a function of the a) photon virtuality squared Q 2 , b) Bjorken 
scaling variable x, c) transverse momentum px and d) pseudorapidity 77. More details in the 
caption of figure @J 
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Table 2: The differential cross-section values as a function of Q 2 , x, pr and i] in the visible 
region defined by 2 < Q 2 < 100 GeV 2 and 0.1 < y < 0.6. The bin ranges, the bin averaged 
cross section values, the statistical and the positive and negative systematic uncertainties are 
listed. 
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Table 3: The differential A cross-section values as a function of Q 2 , x, pt and r]. More details 
in caption of tabled 
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Table 4: The differential K® cross-section values as a function of pif reit and rett in the target 
and current hemispheres of the Breit frame. More details in caption of table [2] 
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- 0.3 


0.55 


0.05 


0.04 


0.04 


0.3 - 1 


0.21 


0.02 


0.03 


0.03 



Table 5: The differential A cross-section values as a function of jxf eit and xz r ett in the target 
and current hemispheres of the Breit frame. More details in caption of table [2] 
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R(A/K S °) 


Q 2 




R(A/K° S ) 


stat. 


syst. (+) 


syst. (-) 


[GeV 2 ] 










2 - 


2.5 


0.406 


0.025 


0.018 


0.019 


2.5 - 


3 


0.390 


0.020 


0.030 


0.030 


3 - 


4 


0.368 


0.019 


0.020 


0.020 


4 - 


5 


0.366 


0.020 


0.018 


0.019 


5 - 


7 


0.347 


0.014 


0.014 


0.014 


7 - 


10 


0.369 


0.016 


0.018 


0.019 


10 - 


15 


0.367 


0.015 


0.013 


0.014 


15 - 


25 


0.360 


0.016 


0.017 


0.017 


25 - 


100 


0.353 


0.012 


0.017 


0.018 


X 


R(A/K°) 


stat. 


syst. (+) 


syst. (-) 


0.00004 - 


0.0001 


0.405 


0.017 


0.024 


0.025 


0.0001 - 


0.0002 


0.390 


0.014 


0.019 


0.020 


0.0002 - 


0.0004 


0.355 


0.011 


0.012 


0.013 


0.0004 - 


0.001 


0.364 


0.012 


0.013 


0.014 


0.001 - 


0.01 


0.329 


0.011 


0.016 


0.017 


Pt 




R(A/iT°) 


stat. 


syst. (+) 


syst. (-) 


[GeV] 










0.5 - 


0.6 


0.24 


0.02 


0.02 


0.02 


0.6 - 


0.7 


0.268 


0.017 


0.009 


0.010 


0.7 - 


0.8 


0.309 


0.020 


0.015 


0.016 


0.8 - 


0.9 


0.334 


0.014 


0.016 


0.017 


0.9 - 


1.1 


0.402 


0.017 


0.015 


0.016 


1.1 - 


1.3 


0.450 


0.018 


0.015 


0.016 


1.3 - 


1.6 


0.442 


0.016 


0.021 


0.022 


1.6 - 


2.2 


0.485 


0.016 


0.019 


0.020 


2.2 - 


3.5 


0.505 


0.027 


0.031 


0.032 


V 


R(A/K°) 


stat. 


syst. (+) 


syst. (-) 


-1.3 - 


-1 


0.331 


0.018 


0.015 


0.015 


-1 - 


-0.75 


0.330 


0.017 


0.018 


0.018 


-0.75 - 


-0.5 


0.350 


0.018 


0.013 


0.014 


-0.5 - 


-0.25 


0.323 


0.016 


0.015 


0.016 


-0.25 - 





0.311 


0.018 


0.017 


0.017 


- 


0.25 


0.337 


0.018 


0.018 


0.018 


0.25 - 


0.5 


0.389 


0.019 


0.017 


0.018 


0.5 - 


0.75 


0.420 


0.023 


0.014 


0.015 


0.75 - 


1 


0.430 


0.020 


0.017 


0.018 


1 - 


1.3 


0.48 


0.02 


0.02 


0.02 



Table 6: The values of the ratio R(A/K"°) of the differential cross-sections for A baryons and 
mesons as a function of Q 2 , x, pr and i]. More details in caption of table |2] 
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R(A/K S °) 


pBreit ( ar g et 

[GeV] 


R(A/iO 


stat. 


syst. (+) 


syst. (-) 


0.5 - 1 


0.329 


0.009 


0.008 


0.009 


1 - 1.25 


0.449 


0.019 


0.015 


0.016 


1.25 - 1.5 


0.443 


0.023 


0.017 


0.018 


1.5 - 2.5 


0.493 


0.019 


0.020 


0.021 


2.5 - 4 


0.495 


0.043 


0.026 


0.027 


pBreit curren i; 

[GeV] 


R(A/iO 


stat. 


syst. (+) 


syst. (-) 


- 0.6 


0.153 


0.015 


0.005 


0.005 


0.6 - 3 


0.185 


0.017 


0.008 


0.008 


x Breit (^j-ggf 


R(A/K°) 


stat. 


syst. (+) 


syst. (-) 


- 0.45 


0.235 


0.014 


0.008 


0.008 


0.45 - 1 


0.277 


0.013 


0.010 


0.011 


1 - 2 


0.297 


0.014 


0.013 


0.013 


2-4 


0.370 


0.016 


0.014 


0.015 


4-7 


0.460 


0.022 


0.017 


0.018 


7-11 


0.46 


0.03 


0.02 


0.03 


11 - 20 


0.61 


0.04 


0.03 


0.04 


x Breit curren t 


R(A/2f°) 


stat. 


syst. (+) 


syst. (-) 


- 0.3 


0.169 


0.015 


0.009 


0.009 


0.3 - 1 


0.174 


0.017 


0.012 


0.012 



Table 7: The values of the ratio R(A/K®) of the differential cross-sections for A baryons and 
mesons as a function of p® reit and x^ rett in the target and current hemispheres of the Breit 
frame. More details in caption of table [2] 
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R(K°/h±) 




Q 2 




R(K°Jh±) 


stat. 


syst. (+) 


syst. (-) 


[GeV 2 ] 










2.0 


— 


2.5 


0.0681 


0.0002 


0.0028 


0.0029 


2.5 




3.0 


0.0645 


0.0001 


0.0026 


0.0027 


3.0 




4.0 


0.0631 


0.0002 


0.0026 


0.0027 


4.0 




5.0 


0.0674 


0.0002 


0.0028 


0.0028 


5.0 




7.0 


0.0670 


0.0002 


0.0027 


0.0028 


7.0 




10.0 


0.0669 


0.0002 


0.0027 


0.0028 


10.0 




15.0 


0.0696 


0.0002 


0.0029 


0.0029 


15.0 




25.0 


0.0717 


0.0002 


0.0029 


0.0030 


25.0 




100.0 


0.0676 


0.0002 


0.0028 


0.0028 


X 




stat. 


syst. (+) 


syst. (-) 


0.00004 




0.0001 


0.0689 


0.0002 


0.0028 


0.0029 


0.0001 




0.0002 


0.0671 


0.0002 


0.0027 


0.0028 


0.0002 




0.0004 


0.0644 


0.0001 


0.0026 


0.0027 


0.0004 




0.001 


0.0671 


0.0002 


0.0027 


0.0028 


0.001 




0.01 


0.0699 


0.0002 


0.0028 


0.0029 


Pt 




R{K»/h±) 


stat. 


syst. (+) 


syst. (-) 


[GeV] 










0.5 




0.6 


0.0499 


0.0001 


0.0020 


0.0021 


0.6 




0.7 


0.0522 


0.0002 


0.0021 


0.0022 


0.7 




0.8 


0.0572 


0.0002 


0.0023 


0.0024 


0.8 




0.9 


0.0633 


0.0002 


0.0026 


0.0027 


0.9 




1.1 


0.0723 


0.0002 


0.0030 


0.0030 


1.1 




1.3 


0.0773 


0.0002 


0.0032 


0.0032 


1.3 




1.6 


0.0872 


0.0003 


0.0036 


0.0037 


1.6 




2.2 


0.0955 


0.0003 


0.0039 


0.0040 


2.2 




3.5 


0.1020 


0.0003 


0.0042 


0.0043 


V 


R(K° s /h±) 


stat. 


syst. (+) 


syst. (-) 


-1.3 




-1 


0.0700 


0.0002 


0.0029 


0.0030 


-1 




-0.75 


0.0696 


0.0002 


0.0029 


0.0029 


-0.75 




-0.5 


0.0656 


0.0002 


0.0027 


0.0028 


-0.5 




-0.25 


0.0635 


0.0002 


0.0026 


0.0027 


-0.25 







0.0654 


0.0002 


0.0027 


0.0027 







0.25 


0.0616 


0.0002 


0.0025 


0.0026 


0.25 




0.5 


0.0646 


0.0002 


0.0027 


0.0027 


0.5 




0.75 


0.0653 


0.0002 


0.0027 


0.0027 


0.75 




1 


0.0670 


0.0002 


0.0027 


0.0028 


1 




1.3 


0.0721 


0.0002 


0.0030 


0.0030 



Table 8: The values of the ratio of the differential production cross sections for mesons and 
charged hadrons as a function of Q 2 , x, px and rj. More details in caption of table[2] 
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